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We perform the quantitative evaluation of the entanglement dynamics in scatter- 
ing events between two insistinguishable electrons interacting via Coulomb potential 
in ID and 2D semiconductor nanostructures. We apply a criterion based on the von 
Neumann entropy and the Schmidt decomposition of the global state vector suitable 
for systems of identical particles. From the time-dependent numerical solution of 
the two-particle wavefunction of the scattering carriers we compute their entangle- 
ment evolution for different spin configurations: two electrons with the same spin, 
with different spin, singlet, and triplet spin state. The procedure allows to evaluate 
the mechanisms that govern entanglement creation and their connection with the 
characteristic physical parameters and initial conditions of the system. The cases in 
which the evolution of entanglement is similar to the one obtained for distinguishable 
particles are discussed. 
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INTRODUCTION 



Quantum entanglement has undoubtly represented one of the most peculiar aspects of 
quantum mechanics as it can be viewed as the furthest departure of quantum world from the 

HQS 

. Furthermore, in the last years, entanglement has been recognized as 



classical one 



the resource for quantum information processing j^, |5(. In this context a growing interest in 
continuous-variables entanglement is arising for various physical phenomena, such as binary 
collisions between unbound or trapped particles, photoionization of atoms in both weak 



2 



and strong fields, and photodissociation of molecules, where the entangled subsystems are 
spatial separated Q £QQ. For bm a ry eeHisio. a pa rti cu,ar attention has been devoted 
to the dynamical creation of entanglement between the two particles, stemming from the 
possibility of using such a system to investigate the foundations of quantum mechanics 
itself: it has been shown that an almost maximal violation of Bell's inequality can occur 
|lo| . The quest for quantum-computing capable devices has also produced great interest 
in entanglement formation in solid state systems where decoherence effects, mainly due 
to carrier-carrier scattering, play a dominant role [ill . . l^ . In particular these effects 
have been studied in the context of electron transport in semiconductors, where the time 
evolution of the entanglement of spatial degrees of freedom has been calculated numerically 
for scattering events between two distinguishable charged particles by means of the von 
Neumann entropy of the density matrix reduced on the spatial degrees of freedom of one of 



the particles 



11 



]. It should be noticed that the distinguishability of the two particles 
greatly simplifies the system modelling and that it is representative of only few of the 
collision events taking place in semiconductors, as an electron-hole scattering. Indeed the 



notion of entanglement for quantum systems composed of distinguishable particles has been 



widely investigated from the theoretical and experimental point of view |8ul4j. while relevan ; 



questions arise for systems of identical particles (both bosons and fermions) 



16 



17 



18, 



19|. Here, in particular, unlike the distinguishable particles case, the Hilbert space of the 



physical states does not have a naturally selected structure of tensor product, because of 
the symmetrization postulate. For example, in the analysis of a scattering event between 
particles, this implies that the system wavefunction can never be factorized in the product 
of two one-particle wavefunctions. This makes difficult to give a straightforward definition 
of entanglement because of the presence of the correlations due to the exchange symmetry. 
Only recently the problem has been deepl y a nal y zed for fermion systems. Various ap- 



proaches have been proposed in literature 



eepi y a nal y zed 101 

0,0,00, E, 



22j, each of one having its own 



advantages and drawbacks, and no consensus has been reached on which one is the most 



suitable |2j|. Here we follow the one introduced by Schliemann (jja,l2J] and recently reexam 



ined and extended by other authors (20 , 



3,0 



since it allows to use a density matrix in the 



spatial coordinates and to compute the von Neumann entropy by means of a generalization 
of the numerical procedure already developed for the case of distinguishable particles. Such 
an approach is based on an analogous of the Schmidt decomposition for state vectors of 
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two fermionic particles: through an unitary transformation the antisymmetric wavefunction 
is expressed into a basis of Slater determinants with a minimum number of non-vanishing 
terms. This number, known as Slater rank, is a criterion to identify whether a system is 
entangled or not, which involves the evaluation of the von Neumann entropy of the one- 
particle reduced density matrix. In this notation the von Neumann entropy becomes not 
only a measure of the uncertainty due to the impossibility of attributing a definite state to 
each particle of the couple but also of the amount of uncertainty deriving from the indis- 
tinguishability of particles. However it is well known that the quantum correlations related 
to exchange simmetry cannot be used to violate Bell's inequality and are not a resource for 
quantum information processing [^J : they cannot therefore be considered as a manifestation 
of genuine quantum entanglement. 

The basic significance of a quantitative evaluation of particle-particle entanglement comes 
also from the link between quantum correlations and decoherence l|. The latter is, in fact, 
a manifestation of the former when one drops the exact quantum description of the whole 
system, considering a part of it as the environment, thus tracing out its quantum coordinates. 
In our case one can think of one of the particles as representing the environment. As a 
consequence the amount of entanglement determined by the proposed simulative approach 
is a direct measure of the decoherence undergone by a carrier as a consequence of a carrier- 
carrier scattering. 

In the present paper we aim at calculating entanglement of both spatial and spin degrees 
of freedom by studying, via a time- dependent numerical analysis, ID and 2D models of 
a scattering event between two electrons in semiconductors. The two particles, explicitly 
considered as indistinguishable, interact through a Coulomb potential. This approach allows 
us to overcome the hypothesis of distinguishabilty of the particles used in previous works 
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ll^ . In the present case quantum entanglement must be analyzed more carefully in 
order to single-out the effects of the additional correlation induced by indistinguishability. 
We also evaluate the role of the spin degrees of freedom. We believe that our results can 
be a valuable guideline for the use of collision events of identical particles in the quantum 
information processing based on solid state systems and for the estimation of the decoherence 
rate induced by carrier-carrier scatterings. 

The paper is organized as follows. In Sec. |H] we describe the method used to quantify 



entanglement in two-fermion systems. In Sec. IHII we investigate the entanglement dynamics 
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for the case of a couple of electrons travelling in two parallel, single-mode, quantum wires 
(ID-system) and in Sec. UVI the same analysis is performed for a scattering event between a 
free propagating and a bound electron in a single-mode two-dimensional electron gas (2D- 
system). Finally in Sec. E)we comment the obtained results and draw conclusions. 



II. ENTANGLEMENT EVALUATION FOR TWO-FERMION SYSTEMS 

Here we shall briefly describe the notion of entanglement between two fermions and the 
approach used for its evaluation, following the theory formulated by Schliemann l^ . 
Any pure two-fermion state can be written as: 

2N 

\®f) = y]wijala\\0) {ai,a]} = (1) 

id 

where Wij are the elements of a complex and antisymmetric (2N x 2N) matrix W, and 2N 
is the number of modes for each single particle. a\(a,i) is the creation(annihilation) operator 
of a fermion in the z-th mode. The normalization condition is Tr[W"'W] = 1. Now we need 
to introduce the fermionic analog of the Schmidt decomposition indicating that exists 
a unitary transformation U such that 

W = UZU T (2) 

where 

Z = diag[Zi ...Zi...Z N ] (3) 
is a block diagonal matrix with bidimensional blocks Z\ of the type: 



Zi 
-Zi 



We note that Z{ may also be zero. This decomposition is unique and we can write 

N 



\* F ) = Y, z ^4i-Ai - 44-i)|o) (4) 



t 



where the fermionic state is given in terms of new creation operators defined from a % 



Y2j u ij a j- m other words one can think of the unitary transformation U as acting on the 
fermionic operator a\. The number of Zj different from zero is the Slater rank and is closely 
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related to the entanglement: the state \^f) is considered non entangled if its Slater rank 
is equal to 1, while it must be considered a genuinly entangled state if its Slater number is 
greater than 1. In fact it can be easily seen from Eq. (J3J that only in the first case \^f) can 
be obtained by antisymmetrization of the product of the single-particle states. 

It has been shown recently that the von Neumman entropy, commonly used in the context 
of non-identical particles, is still a good quantum- correlation measure for identical particles 
(2^, (2^. From the two-particle density matrix p F = \^ F )(^ F \, we can get the one-particle 
reduced density matrix (normalized to unity) [20] 

i=^|^=(^n, (5) 

In terms of the coefficients Zi the normalization condition becomes Y^/iLi \ z i\ 2 = 1/2, while 
the von Neumman entropy e can be computed from: 

N 

e = -Tr[p f In/] = ln2 - 2 ^ \ Zl \ 2 In 2\ Zl \ 2 (6) 

i=i 

where, as stated before, \zi\ 2 are the eingenvalues of the matrix . For identical particles 
the von Neumann entropy reaches its minimum value e = In 2 in the case of non entangled 
states, whose Slater rank is 1, while, for distinguishable particles, the minimum value of the 
entropy, still corresponding to non entangled states, is zero. This behavior can be explained 
if we consider the real meaning of the von Neumman entropy, that must be interpreted as a 
measure of the lack of knowledge about the state of a quantum system. The value In 2 for non 
entangled two-fermion states measures the uncertainty arising from the indistinguishability 
of the particles; in this case the correlations are only related to exchange simmetry and 
cannot be used as a resource for quantum information. On the other side when the von 
Neumann entropy is greater than In 2 there is an additional ignorance about the state of one 
of the two particles, stemming from genuine entanglement. 

III. ID-SYSTEM 
A. The Theoretical Model 



The physical system consists of two electrons running in opposite directions along two 
parallel silicon quantum wires. The dynamics of the two particles is ID and is coupled 
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through Coulomb potential, so that the Hamiltonian of the system takes the form 



where e is the electron charge, e is the silicon dielectric constant, x a and Xb the electron 
coordinates along the wires and d the distance in the perpendicular direction y between the 
two wires. We use the approximation of single parabolic band and do not consider spin-orbit 
interaction. 

The two quantum wires have been considered as ID-systems. Each of the two particles 
has another degree of freedom, besides the spin and space one, indicating the quantum wire 
in which it is moving, i.e. we suppose that each carrier is confined along the y direction in 
a wire (see Fig. [IJ. Moreover each of the two carriers is represented at initial time to by a 





where V(r , r^) can be expressed as 



V(r a ,r b ) 




(8) 
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FIG. 1: (Color online) The ID physical system consisting of two parallel quantum wires, where 
two electrons run in opposite directions; d is the distance between the wires. 
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minimum uncertainty wave packet, so that its space wavefunction takes the form 

. 1 / (x — x ) 2 ., \ , . 

^ to) = exp r~iw" + lk ' x ) (9) 



where a is the mean dispertion in position, k = ^2mE k /h with m the effective mass of the 
carrier and Ek its kinetic energy. 

Taking into account both spin and "wire" degrees of freedom, we can consider differ- 
ent initial state configurations. The first quantum state considered is the one having two 
particles with the same spin (spin up): 

|*) = 4^0)|LR)H^)|RL>W ( 10 ) 



V2 ] 

where both the wavefunctions corresponding to the states and \<p) are of the type defined 
in Eq. (jHJ), and the variances of the wavepackets and the distance between their centers are 
such that the Coulomb energy of the system is negligible at initial time. In Eq. f(10j) . the 
ket |L)(|R)) indicates that the particle is localized in left (right) wire, and ||) indicates spin 
up state. 

In order to compute the entanglement as given by Eq. (j^J) we need to evaluate numerically 
the quantum state of Eq. (jlOj) . To this aim the explicit form of the matrix W has to 
be obtained. This is done in the Appendix „ where we show that W for the same-spin 
configuration takes the form 

/ W A -W s \ 
W s -W A 

\ 0/ 



2V2 



(11) 



where W a and Ws are the antisymmetric and symmetric matrices defined in the Appendix 

D 

The second quantum state considered is the one of two electrons having different spins. It 
is not factorizable in a part containing only space and wire variables and in a part containing 
only spin variables, and reads 

|T) = 4^^)l LR )ltI)-l0^)|RL)|!T)Y (12) 



V2 

It should be noticed that for such a state it is possible to make precise claims about the 
spin of a particle whose position-wire state is known. From the calculations presented in the 



8 



Appendix B we get 





( w A 


-W s 


-W s 


w A \ 


1 


W s 


-W A 


-w A 


W s 




W s 


-w A 


-w A 


W s 




\w A 


-W s 


-W s 


w A ) 



(13) 



The last two states considered still describe two electrons with different spins, but they 
can be factorized in a position-wire coordinate term and in a spin term. We can identify 
the singlet spin state 



|^0}|LR) + i0V)|RL) MTD-ilT) 



(14) 



with 



W* 



/ o 



W s 

\w A 






-w A 

-W s 



-W s W A \ 
-W A W s 








(15) 
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and the triplet spin state 



|~> = -(hA0>|LR> 



|0^)|RL) ITD + I4T) 



(16) 



with 



/ W A -W s 



w s = - 



(17) 



\ 

W s -W A 
-W A W s 
\ -W s W A ) 
The triplet and singlet states have relevant properties that have been taken into account in 
evaluating the von Neumann entropy. In fact, as shown in the Appendix B the von Neumann 
entropy of |S) and |$) can be obtained by adding In 2 to the one of the same-spin state l^). 
This property, due to the fact that such states are factorizable, remains true also during 
the time evolution since the Hamiltonian does not include spin terms. This implies that at 
any time t the triplet and singlet can always be factorized in a position-wire term and a 
spin term. Therefore we can attribute the offset of the entanglement value to the lack of 
knowledge about the spin state, in the sense that one cannot make precise claims about the 
spin state of the particle in a given position-wire state \i/jL) or \<pR) 



9 



B. Numerical results 

In our approach we solve numerically the time-dependent Schrondiger equation for the 
two-particle spatial wavefunction by means of a Crank-Nicholson finite difference scheme. 
Once the spatial wavefunctions are obtained at each time step, we calculate the matrices 
W's related to the different initial spin state configurations and, from these, we obtain the 
one-particle reduced density matrixes p's. This allows us to evaluate the time evolution of 
the entanglement between the spatial and spin degrees of freedom of the two electrons in 
terms of the von Neumann entropy according to the expression given in Eq. (JHJ). To this 
purpose we need to diagonalize the one-particle reduced density matrix at each time step 
and this turns out to be the most demanding part of our algorithm from the point of view 
of the numerical calculation. 

Figs. El and H3 show that, at initial time, when the two particles are far enough and 
Coulomb interaction is negligible, the von Neumann entropy for the states \^f) and |T) is 
In 2 as expected j^,!^ . As indicated by the theoretical criterion described in Sec. ITTlthis 
value can be related to the unavoidable correlations due to exchange symmetry and does 
not imply genuine quantum entanglement. In fact \^f) and |T) are initially non entangled. 
The entanglement increases while the two electrons are approaching each other and finally 
reaches a stationary value once electrons get apart. The final entanglement, which is a 
function of the distance d between the two wires, is reported in the insets of Figs. El and 
01 Following this approach it is possible to obtain a quantitative evaluation of the time of 
entanglement formation in a scattering event: for the studied case such a time is found to 
be around 1 picosecond. Obviously the dynamics of entanglement formation depends on 
the physical parameters of the system, as the initial kinetic energy of the electrons and the 
variance a of their initial wavepacket. This is illustrated in Fig. 0] where e(t) is reported for 
different values of the variance a. Increasing the mean dispersion in position leads to an 
increase in the time at which the stationary value of entanglement is reached, but the value 
itself is pratically not affected. 

We note that the dependence on d of the final values of the entanglement is significantly 
different for the case of the state with respect to the case of the state |T). While for 
this value decreases monotonically as the distance between the two wires is increased, for |T) 
it shows a well defined maximum (around d — 1 nm in Fig. 01 inset). This means that, in the 
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d = 0.50 nm 
d= 1 nm 
d= 1.50 nm 
d = 2 nm 
d = 3 nm 



0.5 1 

Time (ps) 

FIG. 2: (Color online) Entanglement as a function of time for the state desribing two electrons 
with the same spin for different values of the distance d between the two wires. The two electrons 
have the same initial energy = 50 meV and are described by two wavepackets with an initial 
mean dispersion a = 20 nm. The inset shows the stationary values of the entanglement as a 
function of d. 

case of electrons having the same spin, the scattering creates strong quantum correlations 
between the particles which are weakly affected by the wire distance. On the other hand for 
electrons with different spins we find a strong dependence on d of the stationary values of 
the entanglement. This behavior is in good qualitative agreement with the one obtained in 
previous works in the case of distinguishable particles. 

In Fig. we report the time evolution of the entanglement for two values of the distance 
between the wires and for all the considered spin configurations. Intially the von Neumann 
entropy of the singlet |<&) and triplet |H) states has the value 2 In 2. In agreement with the 
theoretical predictions its time evolution results to be the same obtained for the same-spin 
state |$), just shifted by In 2. Thus at the initial time the singlet and triplet states have to be 
considered as entangled according to the notion of entanglement for fermionic particles used 
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FIG. 3: (Color online) Entanglement as a function of time for the state |T) (see Eq. (|12|l) for 
different values of the distance d between the two wires. The two electrons have the same initial 
energy = 50 meV and are described by two wavepackets with an initial mean dispersion a = 
20 nm. The inset shows the stationary values of the entanglement as a function of d. 



in this work. The appearance of the offset above mentioned can be explained considering 
the mathematical form of |$) and |H) as given in Eqs. ()14j) and ()16|) . Their Slater rank is 
greater than 1 since they cannot be put in terms of a single Slater determinant: in fact they 
are not obtained by antisymmetrization of the product of single-particle states. This also 
confirms that the von Neumann entropy, as a measure of entanglement, corresponds to the 
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Slater rank criterion 



IV. 2D-SYSTEM 



A. The theoretical model 



We study in this Section the entanglement formation for the case of an electron propa- 
gating in a 2D system and interacting with another electron bound to a specific site by a 
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FIG. 4: (Color online) Entanglement as a function of time for the same-spin state \^f) at different 
values of the initial mean dispersion a of the spatial wavefunctions. The initial kinetic energy 
is 50 meV for both electrons. Note that the final stationary value is pratically independent upon 
a. 

harmonic potential. The latter can be considered as a simple model of a shallow impurity 
in a semiconductor. Here the physical parameters of GaAs have been used. The dynamics 
of the particle bound in the harmonic potential is coupled to the incoming particle through 
a Coulomb interaction. The Hamiltonian reads 



where e and m are the GaAs dielectric constant and effective mass, respectively, and r is 
the center of the harmonic potential, with energy-level spacing f\w. 

At the initial time to one of the two particles (incoming electron) is represented by a 
minimum uncertainty wavepacket ip(r,to) centered in ri, with k vector pointing towards ro, 
while the other (bound electron) is in the ground state of the 2D harmonic oscillator 
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FIG. 5: (Color online) Entanglement as a function of time: two electrons with the same spin 
| \P) (dashed lines), with different spin |T) (dotted line), singlet |<£>) and triplet |S) (solid line). The 
two wavepackets describing the electrons have both an initial energy = 50 meV and initial mean 
dispersion a = 20 nm. Two different values of the distance d between the two wires are considered: 
d = 0.5 nm in the left panel and d = 2 nm in the right panel. 



The distance |r^ — r | is such that at the initial time t Q the Coulomb energy is negligible. 
Moreover we have taken a = a/ h/2mu) : this allowed us to evaluate the time evolution of 
the two-particle wavefunction in the center-of-mass coordinate system. 

Let us now consider the state describing the system with the two particles having 
the same spin. At the initial time 

i*> = -^^0>-i^>)m>. (20) 

When we discretize each space dimension with a grid of N points we obtain 

N 2 



*> = 7f E (W*)<^) -^(^^(r^lr^lTT) 



(21) 



and applying the same formalism used for the ID system of Section IIHI we get for the 
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2N 2 x 2N 2 matrix W associated to the same-spin state 

1 (w A o\ 

W * = -7k > ( 22 ) 

V2 I I 

where Wa is the antisymmetric A^ 2 x A^ 2 matrix whose elements read 

Wij = ip(r t ) (f>(Tj) - 0(ri) ^(rj), (23) 

obtained by the antisymmetrization of the product of the single-particle wavefunctions ip 
and <fi. The other spin configuration we study in this Section is the triplet spin state 

|H) = i^0)-|0^)V|TI) + |!T)l (24) 



2 



for which 



1 W A . 
Wh = , • (25) 

2 \ -W A 

Also in this case the von Neumann entropy of the triplet |S) is obtained by the one of 
the same-spin state \^f) adding In 2. The singlet state and the state with the electrons 
having different spins and non factorizable into a spin term and a space term have not been 
considered in the 2D case, since they need a much larger computational effort that is beyond 
the scope of the present work. 



B. Numerical results 



The time evolution of the 2D two-particle wavefunction \^f) has been performed in the 
center-of-mass coordinate system, by means of a finite-difference solution of the correspond- 
ing time-dependent Schrodinger equation 30] . 

At initial time the centers of the wavepackets are distant enough so that the Coulomb 
energy is negligible: as a consequence, when the system is described by the state \^/) (the 
two electrons have the same spin), no correlation is present apart from the one due to 
the exchange symmetry. In fact we see from Fig. El that for any value of the energy of 
the incoming electron the initial von Neumman entropy gets its minimum value In 2. As 
the particles get closer their quantum correlation builds up and entanglement reaches a 
stationary value once particles are far enough. Here the final amount of entanglement 
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depends on the initial energy of the carrier: it is higher for higher energies. The time 
evolution of e(t) for the triplet state can be obtained simply shifting by In 2 the curves 
reported in Fig. |H1 

Fig. 0is intended to give a better insight into the role played by the correlations between 
the space degrees of freedom in the entanglement creation, for the |\&) state case. To this 
purpose we have calculated the square modulus of the wavefunction, at two different times, 
keeping the position of one of the two electrons fixed (indicated by the black spot in the 
figures). The first time (left column) corresponds to a condition of minimum entanglement, 
the second (right column) to a condition in which entanglement has already reached its 
stationary values. We note that when entanglement is at its lowest level the result is in- 
dependent on the choice of the fixed position while this is not true once entanglement has 
been created by the Coulomb interaction. 
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FIG. 6: (Color online) The time evolution of entanglement for \^f) (see Eq. (|20)0 in the 2D system 
at three different values of the initial energy of the incoming electron, namely = 10 meV (dashed 
line), 20 meV (dotted line), 30 meV (solid line). The harmonic oscillator energy is frw = 2 meV. 
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FIG. 7: (Color online) The square modulus of the two-particle antisymmetrized wavefunction 
■0(r) 4>(R,) — </>(R) iP(t) for |^) at two different times: t = 40 fs and t = 400 fs. In the two upper 
graphs the position of one of two particles is fixed to R=(100 nm, 100 nm), while in the two lower 
graphs R=(50 nm, 150 nm). 



V. CONCLUSIONS 

In this paper we have studied the dynamics of entanglement formation in systems com- 
posed by two indistinguishable particles. To this aim we applied a theoretical method 
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24j . which permits the quantitative evaluation of entanglement for a two-fermion sys- 
tem. The method is based on a fermionic analog of the Schmidt decomposition of the state 
vector and agrees with von Neumman entropy of the single-particle reduced density matrix, 
which is a good correlation measure and allows one to determine whether the lack of knowl- 
edge is due to a genuine entanglement or only to particles indistinguishabilty. It should 
be noticed that the correlations due to exchange symmetry cannot be used to implement 



quantum information processing or violate Bell's inequality 
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The time-dependent numerical implementation of the theoretical method we developed 
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allows one to investigate the decoherence effects that are consequence of carrier scattering, 
and that can be quantified in terms of entanglement between the scattered particles. This 
approach makes feasible a direct quantification of the time of the entanglement formation. 
The numerical results have been obtained using the effective mass of electrons in semicon- 
ductors (Si for the ID model and GaAs for the 2D model) but they can be considered as 
representative of a general behavior. 

Our simulations show that the final value of entanglement between two electrons un- 
dergoing a scattering event depends only upon some specific physical parameters, as the 
distance between the quantum wires (in the ID case), or the initial kinetic energy of the 
particles, while it is essentially immune to the modifications of other parameters, like the 
variance a of the initial wavepackets. We showed how the entanglement dynamics depends 
on the spin components of the state vector describing the system at the initial time, even if 
the Hamiltonian does not contain spin terms. This suggests the existence of a deep relation 
between spin states and space entanglement, which implies the possibility of detecting quan- 
tum correlations in spin measurements even for the case of wavepackets extending in two 
distant, not overlapping, regions of space. Moreover, in the ID system, the numerical results 
for singlet and triplet spin states shows, at the initial time, quantum correlations related 
not only to exchange simmetry but also to genuine entanglement. They can be attributed 
to the lack of knowledge about the spin state of a particle in a specific wire, in agreement 
with theoretical predictions. 

The results obtained in the case of two electrons havin g di fferent spins are in good quali- 
tative agreement with the case of distinguishable particles |llLll2| apart from the offset term 
deriving from exchange simmetry. This also confirms that the theoretical framework given 



in Refs. |l8 . 
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26] is a suitable criterion for the evaluation of the entanglement in two- 
fermion systems and is equivalent, in the case of distinguishable particles to the commonly 
used criterion based on the Schmidt number or on the von Neumman entropy. 
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APPENDIX: DETERMINATION OF W MATRICES 

Here we illustrate the procedure we apply to the different initial states in order to evaluate 
numerically the von Neumann entropy in ID physical systems with each electron localized in 
the left (|L)) or right (|R)) wire, as explained in Sec lIIII Let us begin with the state \^f) defined 
in Eq. (|10|) where both particles have spin up but different position-wire wavefunctions 
described by or |0)|R). We discretize the x coordinate of the electrons along the two 

wires with a N points grid: the dimension of each single-particle Hilbert space is 2N, taking 
into account the spin degrees of freedom. In this way we can consider a fermionic position 
operator at(xi) that creates a particle in the point xf. a^(x/)|0) = \xi). 

We can then express the state |\]/) as 

I*) = ^e(^)^wi lr >-^)^)i rl >) 

x |x/x m )|tT). (A.l) 

We observe that in Eq. (jA.l|) the state is not antisymmetric with respect to the exchange 
of the I and m indeces, and, in this form, it cannot be used to evaluate the entanglement 
according to the theory developed in Sec. HU This difficulty can be overcome by applying a 
unitary transformation to the wire degrees of freedom: 

\V X ) = 1/V2(|L} + |R) j 

> < (A.2) 

\V 2 ) = 1/V2ML) - |R)J. 

Inserting the closure relation in Eq. (jA.l|) . we get 

1 N 2 ( \ 
1^ = 7f5Z Yl ( ^(^)0O£m)7L(s)7RW - <i>(x l )il)(x m )'y R (s)'y L (t) J 

* l,m s,t=l ^ ' 

x |x^ m )|TO)|TT) (A.3) 

where we introduced 7l(1) = 7l(2) = 7r(1) = ^ and 7r(2) = — 

Now we can express the state |^) in Eq. ()A.3|) in terms of fermionic operators fj defined 

as: 

|xi)|1/i)|t) for l<i<N 



//|0) 



|^_7v)|^ 2 )|T) for N + l<i<2N 
\xi-m)\Vi)\l) for 2iV + l<i<3iV 
\xi-w)\V 2 )\l) for 3N + l<i<4N 
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and we get 

47V 

*> X>^4l°> (A.4) 

where the 4iV x 4A" matrix of the coefficients Wij takes the form of Eq. (jlljl . where Wa is 
an antisymmetric N x N matrix whose elements read 

WA lm = tp(xi)4>(x m ) - (j){xi)ip(x m ) (A. 5) 

and Ws is its symmetric conterpart. They have been obtained by antisymmetrizing or 
symmetrizing the product of the single-particle spatial wavefunctions and ip calculated in 
the grid points I and m. 

Now the whole matrix Wq, is antisymmetric and the state in Eq. (|A.4j) is expressed 
in a form that allows us to apply the Schliemann theory The normalization condition 
TrfVy^W 7 ^,] = 1 is satisfied. Moreover during the time evolution Wq, mantains this form 
since the Hamiltonian does not contain spin terms and the coefficients corresponding to the 
modes where particles have spin down will always be zero. 

For the initial state |T) of Eq. (|12jl . with two electrons having different spins and not 
factorizable in a part containing only space and "wire" variables and in a part containing 
only spin variables, we again discetize the x coordinate into a N points grid and consider the 
change of basis given in Eq. (jA.2|) . In this case we have to apply a unitary transformation 
also for the spin variables: from ||) and ||) to \Si) and \S 2 ) 

|^> = i/V5(|T> + ID) 
\s 2 ) = i/V2(|T> - I!)). 

Now we can express |T) as 

l,m s,t,u,v=l ^ 

0(x/)^(x m )7 R (s)7 L (t)/3 i (M)/?t( t ') ) \xix m )\V s V t )\S u S v ) (A.7) 



(A.6) 



where = ft(2) = &(1) = ^ and ft(2) = 

If we define the fermionic creation operators g\ as: 



\xi) |VL)|Si) for l<i<N 

|^_/v}|y 2 )|5'i) for N + l<i<2N 

Xi-2N)\Vi)\S 2 ) for 2N + l<i<3N 

{ \xi- m )\V 2 )\S 2 ) for 3iV + l<«<4iV 
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the state |T) in Eq. (JT2J) can be expressed as 

47V 

|T> = ^^J|0>. (A.8) 

The AN x AN matrix W r reads as in Eq. (|T3j). 

Finally we are able to apply to the singlet spin state of Eq. (fTlj) and to triplet spin state 
of Eq. (jTfij) the same procedure used above for |T), and we find that |$) is given in an 
expression similar to Eq. (|A.8|) . for which the matrix W takes the form of Eq. (j!5jl . while 
for the state |S) W takes the form of Eq. (fTT)|) . For the sake of brevity we do not report 
explicitely those calculations. 

Now we shall show that the von Neumann entropy £= of the triplet |H) and the singlet 
|$) states (both representing the electrons with opposite spin) can easily be related to the 
entropy of the same-spin state Firstly we note that |H) and |<£>) have the same 

entropy: this depends on the fact that the one-particle reduced density matrix W^W results 
to be the same in the two cases. Secondly we observe that the matrix W3 can be viewed as 
a block diagonal matrix with 2 blocks of dimension 2N. Therefore for the eingenvalues \zi\ 2 
of WiW^ holds the relation 

\4? = kf+2ivl 2 - (A.9) 

In this case the e= can be written as: 

2N 

£ E = \n2-2^2\zf\ 2 \n2\zf\ 2 . (A.10) 

The matrix has the same upper left block as and is zero elsewhere. This implies that 
the first 2N eingenvalues of the matrix W^W^ are equal (apart from the constant factor \ 
coming from the normalization condition) to the first 2A^ eingenvalues of the matrix WjWs. 
Using (JA.9)) and (JA.10)) we can relate the von Neumann entropy for |S) to the one for |^) 

- 1 1 

£s = ln2-4^-|zf | 2 ln2-|zf| 2 = e* + ln2 (A.ll) 

i=i 

where we have taken into account the normalization condition Yli=i \ z f\ 2 = |> recalling 
zf = for i > 2N. Eq. (jA.lljl implies that the von Neumann entropy of the triplet spin 
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state can be simply obtained from the one of same-spin state by adding In 2. 
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